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Про моделі промислового 
циркулярного нагнітання 
цементного розчину в сухий грунт 


Представлено моделі числового розрахунку еволюції в часі та просторі фронту промислового циркулярного 
нагнітання цементного розчину в сухий грунт. Першою метою даної роботи є демонстрація того, що 
криволінійна сітка, на якій виконується цей розрахунок, на деяких часових шарах може мати хаотичне 
розміщення вузлів у просторі. Іншою метою є перевірка того, що цей факт не призводить до неузгодженості 
результатів числових розрахунків згідно з представленими моделями. 


Вступ 


Перед будівництвом гребель та інших гідротехнічних споруд часто виникає не- 
обхідність у покращенні механічної міцності їх основ та у створенні протифільтраційних 
завіс. Цього можна досягти нагнітанням через інжектор у грунтовий масив в'яжучого роз- 
чину. Відомо два варіанти нагнітання: циркулярний та затискний. При циркулярному 
способі тиск нагнітання є постійним, а при затискному інфільтрат поступає в інжектор 
при постійній витраті. 

Силікатні гелі, які раніше використовувались як інфільтрат, були хімічно-нестабіль- 
ними і не задовольняли екологічним нормам. Тепер вони замінені цементними розчи- 
нами. Останні складаються з частинок. У випадку, коли їх розмір є достатньо малим для 
того, щоб рухатись у грунті, вони є відповідальними за явище, яке в літературі отримало 
назву глибинна фільтрація або фільтрація безодні. Розрізняють два механізми такої фільт- 
рації. Згідно з (1| перший відповідає затисненню частинок у вузьких проходах між по- 
рами трунту, в той час як другий механізм відображає поступове утворення осаду з цих 
частинок на стінках пор та у проходах між порами, що має наслідком поступове змен- 
шення розміру пор. У моделі нагнітання цементного розчину в насичене водою пористе 
середовище |2| використовується перший механізм для опису фільтрації безодні, а поста- 
новка задачі та значення вхідних параметрів відповідають умовам, наближеним до проми- 
слового використання цієї технології. Результати розрахунків згідно з останньою модел- 
лю свідчать про те, що за цих умов, глибинна фільтрація не є домінуючим явищем. 

У роботі |3| показано, що при виконанні умови 


/ чо» ар, И) 


де тут і надалі Й - швидкість частинок рідкої фази, а а; - коефіцієнт повздовжньої 
дисперсії, описуючи нагнітання цементного розчину в насичений водою трунт, можна 
знехтувати особливостями поширення розчину в пористому середовищі і моделювати 
даний процес задачею з вільною рухомою межею. Фільтрація у грунті цементного роз- 
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чину є прикладом повзучого руху (4). При такому русі силами інерції можна знехтувати 
порівняно з силами тертя. Оскільки цементний розчин є структурованим 1 комбінація 
сил взаємодії між твердими частинками, між твердими частинками та частинками 
рідини робить його більш однорідним зі зростанням концентрації (51, то часткове запов- 
нення ним пор грунту при його нагнітанні в сухий грунт можна моделювати як його 
розчинення в фіктивній невагомій рідині нульової в'язкості, яка насичувала грунт 
перед нагнітанням. Явище механічної дисперсії при нагнітанні цементного розчину в 
насичений рідиною грунт (а / 2» уоде бо коефіцієнт дифузії) пояснюється тим, 
що швидкість частинок речовини в порі грунту плавно змінюється від нуля на поверхні 
до деякого максимального значення в деякій внутрішній точці пори 1 тим, що пори у 
грунті утворюють хаотичну систему з'єднаних між собою тунелів |6|. Якщо порівню- 
вати при інших рівних умовах нагнітання цементного розчину в насичений водою грунт 
з нагнітанням цементного розчину в сухий грунт, то в останньому випадку відношення 
кількості частинок цементу, які знаходяться одночасно в зоні переходу від грунту з 
максимальною концентрацією цементу в порах до грунту з нульовим значенням цієї 
концентрації і в контакті з поверхнею скелету грунту, до повної кількості частинок 
цементу, які знаходяться в зоні цього переходу, є більшим, ніж в першому випадку, 1 
тому при моделюванні сухого грунту в останньому випадку грунтом, насиченим фіктивною 
невагомою рідиною нульової в'язкості, значення а; є меншим, ніж відповідне значення в 
першому випадку |3|. Тому при перевірці виконання умови (1) для випадку нагнітання це- 
ментного розчину в сухий грунт в її праву частину замість а; можна підставляти значення 
для насиченого водою грунту. В роботі |3| показано, що при промисловому нагнітанні 
цементного розчину в сухий грунт умова (1) виконується. Моделі промислового цирку- 
лярного нагнітання цементного розчину в сухий грунт |З3|, (71, І9| формулюються у криво- 
лінійному чотирикутнику |З9| із змінною в часі і наперед не відомою формою, і тому для їх 
розв'язання використовується різницевий метод з використанням числових конформних 
відображень |3), 1101. У цьому методі використовується алгоритм генерації на кожному 
часову шарі криволінійної сітки, яка покриває цей чотирикутник. 

Першою метою даної роботи є демонстрація того, що на деяких часових шарах 
ця сітка може мати хаотичне розміщення вузлів у просторі. Іншою метою є перевірка 
того, що цей факт не призводить до неузгодженості результатів розрахунків згідно з 
моделями |3), (71, 181. 


Моделі промислового циркулярного нагнітання 
цементного розчину в сухий грунт 


У роботах |31, (71, І3| ми розглядаємо чотири постановки задачі. Для постановки 
Хо 1 |7) та постановки Хо 3 |3|, |8| вважаємо, що під фундаментом інжектора зроблено 
довгу траншею ширини 2 - 7; та глибини Ло, в яку при постійному п'єзометричному на- 


порі Но нагнітається (рис. 1) в'яжучий інфільтрат, який є нестисливою ньютонівською 


рідиною. У постановках Хо 2 |7| та Хо 4 |3|, |8| замість траншеї маємо свердловину радіуса 
го та глибини /д. У всіх цих задачах необхідно знайти еволюцію в часі та просторі фрон- 
ту нагнітання. У перших двох постановках скелет грунту вважається абсолютно твердим, 
а у третій та четвертій постановках він є пружним. 

Ми припускаємо (гіпотеза Хо1), що при протяжностях процесу та при тисках в 
інжекторі, які відповідають промисловому використанню цієї технології, початкова 
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форма рухомої межі Г; (рис. 1) не впливає на її кінцеву форму, якщо перша несуттєво 
відрізняється від форми інжектора 1 що у випадках третьої та четвертої постановок 
початкове значення функції напору можна вибрати довільним. Таким чином початкове 
положення рухомої межі Г; вибираємо так, щоб воно несуттєво відрізнялось від Г3. 


Рисунок І - Криволінійні чотирикутники С(г) та С(г) 
Вважаючи фронт нагнітання інфільтрату у грунт, насичений фіктивною невагомою 
рідиною нульової в'язкості, лінією розриву параметрів фільтраційного потоку, легко от- 
римати закон за яким рухаються точки фронту нагнітання |7| 


бентбр, о РО, ой (т гейоо - 


де г (г) - радіус-вектор точки фронту нагнітання в момент часу /, К - коефіцієнт фільт- 
рації, т -- пористість грунту, г - крок по часу, пд - зовнішня відносно до області 
закріплення нормаль до Г, в момент часу / в точці з радіус-вектором 7 (г) : 

У випадку задачі, яка відповідає постановці Хо 1, в це рівняння замість лСа змі 
потрібно підставляти розв'язок крайової задачі |7|, яка формулюється у криволінійному 
чотирикутнику С(г), який на рис. 1 обмежений кривими Г|, Г», Гз та Г,. Для знахо- 
дження розв'язку цієї крайової задачі застосовуємо різницевий метод із використанням 
числових конформних відображень. Цей метод пропонує алгоритм знаходження в кож- 
ний момент часу / - 4, 2 0 такої заміни змінних ху 2 х|(2,1,4), хо з хо (2.п,-), яка 
конформно відображає область С(и ) на параметричний прямокутник С"(п) зі сторона- 
ми Й М (п) та 1, де М (п) - модуль прямокутника (рис. 2 а). Крім того, до розв'язання 
вищезгаданої крайової задачі в момент часу / - | 2 0 застосовуємо метод становлення, 
тобто розглядаємо таку початково-крайову задачу 


ди/дї « ди(2,пь.,ї) є д'и/ да --и/дт, 0) 


и 


пУНоги 


ця зобУМ лм), ди/диуу 2 0. ди/дпт, 2 05 Щепц0)- от), (4) 
де ф - довільна функція змінних 2,7 та (, Гі, Го», Г; та Г; є образами відповідно 
Г, Го, Гі та Гд при вищезгаданому конформному відображенні. Розв'язок задачі (3), 
(4) володіє тією властивістю, що ще тк, 1) зи (2, ц) при 7 -з ко, де й (2,п, й ) - 
-- (х, (2. 7,1 ) 1 ра ДЯ ) у) ,а ла хо, й ) - шукана функція. Знайшовши таким чином 
посуха), де (аоаїе Сс(а), згідно зі співвідношенням (2) знаходимо положення 
фронту нагнітання Г, на момент часу І; - 7, після чого аналогічно знаходимо 


ЛЕО чт) при ехо є си чт), і так діємо до кінцевого моменту часу. 
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Математична модель задачі, яка відповідає постановці Мо 2, отримується з вищеопи- 
саної моделі заміною рівняння (3) на таке рівняння 


ди/дї « д'ші/дЕ? чн (дх» / дп ди/ де - дх» /дЕ - ди/дп)/х (Є трі) д?и/дп?. (5) 


У випадку задачі, яка відповідає постановці Хо 3, в рівняння (2) замість лСа ,Х» .1) 
потрібно підставляти розв'язок крайової задачі (31, яка формулюється у криволінійному 
чотирикутнику С(г). Для знаходження розв'язку цієї крайової задачі так само, як в 
попередніх двох моделях, застосовуємо різницевий метод з використанням числових 
конформних відображень. Однак в цьому випадку більш зручним є нерухомий квадрат 
О із стороною, рівною І ( рис. 2 б). До нього від вищезгаданого прямокутника СТ(І) 
можна перейти за допомогою відображення 2 «2/М (г ) ,пе7,1-1.У роботі |3| по- 
казано, що знаходження функції (ху, хо, ) при (хі, х» ) є С(л ) та ц 2 с.» 0 зводиться 
до розв'язання таких двох одновимірних крайових задач 


ди/дї « Ко(Єп,ц - т)-д'и/дт? ««КЕП,ц -т)-ди/дТ, О« 2,7 «1, | (б) 


Щек -т)е Ер - т), ди/д3(5.0,1)- 0, ди/д5у (2.,7)- 0, (т 
ди/ді а КЕ тп - т) дм дЕ? КОТ, -т)-ди/д2, «2, раї, || (8) 
"ИЙ -г)ецє й), 0,7) є Но, Ії) 2 -хо (117,1). (9) 


Явні вигляди функцій К; РИ -т), де із14, наведені в роботі |31. У рів- 
няннях (6) - (9) п - т 41 «ц,ав рівняннях (7), (9) 042,7 41. У роботі | 3) доведено 
таку теорему. 

Теорема 1. Нехай існують такі функції ц(2,,77,ї) та м(2,7,1) та існує таке б » 0, що 
ці функції є двічі неперервно-диференційованими на множині 0; «(-6,1- 6)х(- 6,16) 
х(д -б- т кб) і є розв'язками задач (6), (7) та (8), (9) відповідно. Нехай існує така 
функція й (2 ОЇ ) , яКа є двічі неперервно-диференційованою на тій самій множині і яка 


є такою, що (туя но (277,0) хо (2,0) при (2т.ї)є О 5 , де (хі, хо, і) - шу- 


неперервними на о,. Тоді можна стверджувати, що виконується наступне: у/ - гу ог? ), 
при і-ї та 044, 141, у-й при їз, 2-0, 0471 та 2-1, 04741 


2 


ди/дії « Ос?) приїоп, їе0,047х1та ло, 42 


п г; 


- 
| зі 
о жи яю 


а) 5) 


Рисунок 2 - Параметричний прямокутник СТ(І) та одиничний квадрат 0 
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Знайшовши таким чином охо, де (ці,ха)є С(ца), згідно зі співвідношен- 
ням (2) знаходимо положення фронту нагнітання Г, на момент часу / -т, після 
чого аналогічно знаходимо Л(ху,х»,-ют) при (х,х») є С(г - т), і так діємо до 


кінцевого моменту часу. 
Математична модель задачі, яка відповідає постановці Хо 4, відрізняється від 


попередньої моделі тільки видом функцій К; (є ок , де і - 3,4, явний вигляд яких 
представлений в роботі |31. 

У кожній із сформульованих вище моделей при розрахунку еволюції фронту 
нагнітання на кожному часовому шарі (: апа 0) область закріплення СІ) покри- 


вається рідкою криволінійною сіткою, наслідком чого є те, що похибка апроксимації 
є суттєвою. У кожній з цих моделей також робляться такі припущення. По-перше, 
ми припускаємо (гіпотеза Мо 2), що при йо »» по в межах цієї похибки фронт 


нагнітання можна описувати двічі неперервно-диференційованою функцією змінної 
х» (рис. 1). По-друге, ми припускаємо (гіпотеза Мо 3), що в межах цієї похибки у 
криволінійного чотирикутника си ) , сторонами якого є криві Г;, де і -1,4, всі кути 


при вершинах є прямими. По-третє, внаслідок того, що кількість вузлів є незначною, 
ми припускаємо (гіпотеза Хо4), що хаотичність на деяких часових шарах розміщення 
у просторі вузлів сітки не є причиною суттєвого спотворення кінцевого положення 
фронту нагнітання. У роботі |3| коректність гіпотез Мо 2 - Хо 4 перевірено для 
постановок Хо» 3 та Хо 4 наступним числовим експериментом. Оскільки інформація 
про те, що фронт нагнітання змінює своє положення, поширюється за допомогою 
звукових хвиль, які внаслідок тертя між грунтом та інфільтратом швидко затухають, 
то можна вибрати криву в межах закріпленої області і момент часу /0, починаючи з 


якого значення функції напору в точках цієї кривої не залежать від часу. Тому можна 
модифікувати моделі для третьої та четвертої постановок, представлені в роботі |З, 
використовуючи припущення про те, що в кожному випадку функція напору 


залежить від часу тільки в області С(і ), яка на рис. 1 обмежена кривими Г,, Г,, 
Г,, Г..1 тому, починаючи з деякого моменту часу 1/0, у рівняння (2) замість функції 
напору підставляємо функцію, яка при 4 2 Іо - 7 в кожному випадку є розв'язком 
крайової задачі |3|, яка формулюється у криволінійному чотирикутнику С(ц). При 
по2 іо ж т модифіковані моделі отримуються з відповідних не модифікованих замі- 
ною відображення одиничного квадрата О на криволінійний чотирикутник Сг) на 
відображення цього квадрата на криволінійний чотирикутник С(и ) та заміною другої 
з умов (9) на умову и/0,77,/)- (0,715), де ке ло)е ЛЕ (2 ук), (а ет 
вбоухомго) - функція напору, отримана в рамках відповідної не модифікованої 
моделі, а ху з Хі (є 7 іо), хр зх» (є з 15) - вищезгадане відображення, яке відобра- 
жає одиничний квадрат О на криволінійний чотирикутник Сід ). В обох цих модифі- 


кованих моделях положення Г3 та момент часу і 0 підбираються експериментально так, 


щоб кінцеве положення фронту нагнітання відрізнялось від кінцевого положення, отрима- 
ного в рамках відповідної моделі з роботи |3|, на як можна меншу величину і щоб по- 
хибка апроксимації розрахунку цього положення в рамках відповідної модифікованої 
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моделі була як можна меншою. При цьому виявилось |З |, що як для постановки ХО 3, так і 
для постановки Хо 4 похибка чисельного розрахунку кінцевого положення фронту 
нагнітання в рамках модифікованої моделі є меншою, ніж в рамках відповідної моделі з 
роботи |8|. Цей факт підтверджує коректність гіпотез Хо 2 - Мо 4 для моделей, які 
відповідають постановкам ХО 3 та Хо 4. 

Вектор нормалі пд до кривої Г;, який присутній у співвідношенні (2), обчислюєть- 
ся так само, як в роботі |8|. Для знаходження числових розв'язків задач (3), (4) та (5), (4) 
використовуємо явні різницеві схеми, представлені в роботі |7), а для числового розв'язан- 
ня пари крайових задач (6) - (9) використовуємо різницеву схему, представлену в роботі (81. 


Алгоритм числового відображення криволінійного 
чотирикутника на параметричний прямокутник 


Застосовуючи різницевий метод з використанням числових конформних відобра- 
жень до розв'язання крайової задачі для диференційного рівняння в часткових похідних, 


сформульованої у криволінійному чотирикутнику с (г) зі сторонами ві (1), де із-14,а 
г - час, роблять заміну змінних 
рені), | заежієньї), (1-Ї» 10) 
яка в кожний момент часу 7 конформно відображає параметричний прямокутник сг) 
(рис. 2 а) на криволінійний чотирикутник с"() (рис. 3). Якщо прийняти (рис. 2 а), що 
довжина вертикальної сторони сг) рівна 1, а довжина Його горизонтальної сторони рів- 
на 1/ М(Г), де М(7) - модуль С" (7), який є інваріантом конформного відображення (9), 
то при знаходженні заміни змінних (10) зручно зробити заміну змінних 
г-с/МІ), | пет, Тзі, (1) 
яка в кожний момент часу відображає одиничний квадрат О (рис. 2 б) на пара- 


метричний прямокутник СІ ). Функції х| (21), х» (27,1), записані в нових змін- 
них, шукаються як розв'язок такої задачі 


м? дж 983 зда дп «0, м? да | дЕ? зда дп «б, (12) 
Фа (Є /мІЄ) 7.8) ж (Є /М(8) т.0)г)-0, (13) 
М дж /дЕ з дха | дП або || дх /дпП «-М дхь /д2. 14) 


Гі 


г/9 


Рисунок 3 - Криволінійний чотирикутник с 1) 
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У рівняннях (12) - (14) 
м з Їжу /дт)? (джу /2п? м аБат / Ще, / 22)? (дез /дЕР лават, (5) 
о о 


де / з ду /дЕ -дх» /д1 - дж /дт -бх» /дЕ . У рівняннях (13), (14) (2, 17)єГ,, а Г; є 


прообразом о б З деіз14, при суперпозиції відображень (11), (10). У рівнянні (13) 
Ф,,деі- 14, є таким, що Ф ог) з- 0 є рівнянням сторони Пи (г) криволінійного 


чотирикутника с (г) в момент часу Ї. 


У роботі (10) представлено алгоритм знаходження числового розв'язку дискре- 
тизації задачі (12) - (14). Але в загальному випадку система рівнянь, яка є дискрет- 
ним аналогом рівнянь (12) - (14) і в якій кількість рівнянь рівна кількості невідомих, 
має необмежену кількість розв'язків. Дійсно, якщо ми покриємо одиничний квадрат 
рівномірною сіткою 3х3, то одна з можливих дискретизацій рівнянь (12), першого з 


рівнянь (14) при (2, 17)єГ,, де і - 12 та другого з рівнянь (14) при (2, 17)є Г,, де 


і «34 ,єтакою 


2 
4М ЕРА 2 вза) 4 бю -2х1 зхуа)в 92 (16) 
2 
ФоМ: бої - 2Х241 у хогі)н 4 еяо - 2211 зхола)є 0, (17) 
Е- Мо -Хід0 ) з 3х210 7 4Х211 1 Х212» Маг -Х102 ) з 3212 З 9Х291119Х210» (18) 
х102 7 2100 7 мОЗхаді -4Х211 1 Х221 ), -х122 1321207 міЗхазі -4Х911 1 Х201 . 19) 


Справедливою є така теорема. 

Теорема 2. Система рівнянь (17) - (19) є лінійно залежною. Доведення. Якщо праву 
і ліву частини кожного з рівнянь (19) помножити на М , і потім до суми цих рівнянь 
додати рівняння (18), то отримається рівняння (17). Теорему доведено. Таким чином в 
результаті виконання вищезгаданого алгоритму (10) можна отримати зайвий розв'язок, 
який не відповідає конформному відображенню. Його ознакою може бути хаотичне роз- 
міщення вузлів криволінійної сітки. 


Результати числових експериментів 


У даній роботі ми розрізняємо два типи вхідних параметрів. До першого типу 
відносимо ті вхідні параметри, які визначаються методом обчислення, а до параметрів 
другого типу відносимо решту вхідних параметрів. Значення параметрів обох типів для 
першої та другої постановок взяті з роботи |7|, а для третьої та четвертої з роботи |31. 
При цьому значення параметрів другого типу для випадку першої постановки є такими 
самими, як для випадку третьої постановки. У свою чергу, значення параметрів другого 
типу для випадку другої постановки є такими самими, як для випадку четвертої. Для всіх 
постановок значення параметрів другого типу відповідають режиму промислового ви- 
користання даної технології. У роботі |3| кінцеві положення фронтів нагнітання для най- 
більш твердого грунту та найбільш м'якого грунту розраховані методом |З) у випадках 
третьої постановки відповідно з похибками 813 7.53 Уо та 213 - 5,090, при цьому міра 


різниці |3| між ними є такоюді з - 1,2 о, а у випадку четвертої постановки відпо- 
відно з похибками є) 4 - 7,30 та Єд - 7,0 96, при цьому міра різниці між ними є та- 


кою 014 - 2,090, У свою чергу, кінцеві положення фронтів нагнітання для найбільш 
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твердого грунту та найбільш м'якого грунту розраховані модифікованим методом |З3| у 
випадках третьої постановки відповідно з похибками є» 3 - 3,7 0 та Є» 3 - 3,690, при 


цьому міра різниці |3| між ними є такою д» з - 1,590, а у випадку четвертої постановки 


відповідно з похибками 624 З 4,690 та 224 - 4,290, при цьому міра різниці між ними 


є такою дод - 2,2 9. Оскільки біз «в бі) ,де і-12, 7 - 3,4, то розрахунки для 


пружного скелету грунту можна порівнювати з розрахунками для абсолютно твердо- 
го скелету грунту |111. 

У роботі |7| кінцеве положення фронту нагнітання для абсолютно твердого скелету 
грунту розраховане у випадку першої постановки з похибкою є; - 2,5 70, а у випадку 
другої -- з похибкою є» з 5,3 Ус. Міри різниці між кінцевим положенням фронту нагні- 
тання в абсолютно твердий скелет грунту, отриманим в рамках першої постановки, та 
кінцевими положеннями, отриманими в рамках третьої постановки методом |8| при 
найбільш твердому грунті, при найбільш м'якому грунті, модифікованим методом |3| при 
найбільш твердому грунті, при найбільш м'якому грунті, відповідно є такими у/ з - 3,9 7, 
б 3 7. 3,5 90, моз 2 6190, У)2 3 - 5,190. Міри різниці між кінцевим положенням фронту 
нагнітання в абсолютно твердий скелет грунту, отриманим в рамках другої постановки та 
кінцевими положеннями, отриманими в рамках четвертої постановки методом |З8| при 
найбільш твердому грунті, при найбільш м'якому грунті, модифікованим методом |3| при 
найбільш твердому грунті, при найбільш м'якому грунті, відповідно є такими 


ул д 7 12190, Ул д 2 11690, м/о д - 12,69, Мі» д -- 10,59, Той факт, що му; ; та М); ; Не Є 


суттєво більшими, ніж відповідно є; ; 5 Є; 2 Та Є; КЕ; 2» де із12, а )-34, 
підтверджує коректність гіпотез Хо 2 -- Хо 4 для всіх представлених у даній роботі моделей 


П Щ. Коректність гіпотези Хо» 1 для моделей, які відповідають першим двом постановкам 
та останнім двом, перевірено відповідно в роботах |7| та |31. 


Висновок 


У представлених моделях на кожному часовому шарі просторові координати вуз- 
лів криволінійної сітки, на якій виконується розрахунок еволюції в часі та просторі фрон- 
ту нагнітання, знаходяться з розв'язку системи алгебраїчних рівнянь. У даній роботі 
вперше показано, що ця система в загальному випадку має необмежену кількість розв'яз- 
ків. У даній роботі також вперше чисельно перевірено коректність гіпотез Хо 2 -- Мо 4 для 
тих представлених моделей, в яких скелет грунту є абсолютно твердим. 

У моделях |3), (71, І8| у будь-який момент часу фронт нагнітання інтерполювався 
одним і тим самим двічі неперервно-диференційованим локальним сплайном. Вплив не- 
визначеності у способі інтерполяції фронту нагнітання на похибки розрахунків згідно 
них буде досліджуватись. Крім того, результати цих розрахунків будуть порівнюватись 
з результатами польових досліджень для того, щоб обгрунтувати нехтування фільтра- 
цією безодні при моделюванні промислового нагнітання у грунт цементного розчину. 
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о 


Н.Б. Демчук 

О моделях промьшпленного циркулярного нагнетания цементного раствора в сухой грунт 
Представленьт модели численного расчета зволюции во времени и пространстве фронта промьшиленного 
циркулярного нагнетания цементного раствора в сухой грунт. Первой целью данной работьт есть демонстра- 
ция того, что криволинейная сетка, на которой вьшолняется зтот расчет, на некоторьтх временньїх слоях 
может иметь хаотическое расположение узлов в пространстве. Другой целью єсть проверка того, что зтот 
факт не приводит к противоречивости результатов расчетов согласно представленньтх моделей. 


М.В. Фетспик 

Арбоші Веа! 5саїс Модеїз ої а Сетепі Стоиі Птегіесбіоп аї а Соп5їапі Рге5биге іп а Дгу 5оїії 

Моавіз ої а питегіса! саїсиіайоп ої а йте апа 5расе суоішіїоп ої а Попі оГа геаї зсаїе іпіесбїоп оїа 
сетепі ягоці іп а дгу 501 аге ргезепіса. ТБе Нг5і айт ої із ууогк 15 їо 5пому аг а сигуіййпеаг єгід їБі8 
саїсиїаноп 15 регіогтеа оп сап Баме сраойс дї5розійоп ої їз роїпі8 іп врасе оп 50пе те Іауег5. Апоіфег 
айт 15 їр сбесК їбає Ші8 Касі дое5 пої гезиії іп сопННсі атопе гезиіїв ої саїсціайопя ассогдаїпе їо ргезепіей 
тодсеіз. 


Стаття надійшла до редакції 24.02.2011. 
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